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Stability of the Couette-Poiseuille flow by 
the Reynolds-Orr energy equation 

F. LAM 

The normal-mode analysis of the Reynolds-Orr energy equation governing the stability of 
viscous motion for general three-dimensional disturbances has been revisited. The energy 
equation has been solved as an unconstrained minimization problem for the Couette- 
Poiseuille flow. When the flow is subjected to two-dimensional streamwise disturbances, 
an exact solution for the vorticity can be obtained in the case of plane Couette flow. In ad- 
dition, it is found that there exists a minimum Reynolds number (R m in = 42.32) among 
the Couette-Poiseuille profiles. This Reynolds number is lower than the value by Orr 
(44.3) for plane Couette flow. Similarly, for two-dimensional spanwise disturbances, the 
lowest Reynolds number (R m in = 19.88) has been found. For fully three-dimensional dis- 
turbances, it is shown that the minimum Reynolds number is in general smaller than the 
corresponding two-dimensional counterpart for all the Couette-Poiseuille profiles except 
plane Couette flow. In particular, the minimum Reynolds number for plane Poiseuille 
flow is 47.2 and it is associated with a three dimensional disturbance whose spanwise 
wavenumber is about three times of the streamwise wavenumber. The present results 
show that the Squire's theorem does not hold for the basic flows other than the linear 
profile. 



1. Introduction 

Consider a basic flow of an incompressible viscous fluid with velocity vector U, and 
let v! denote the velocity vector of the perturbation. The combined flow field U + vl is 
still solenoidal and its motion must be go verned the Navie r-Sto kes equatio ns and subject 
to the same boundary conditions as U. iRevnoldsl f|l895h and [53 (I1907L pp. 122-138) 



showed that the rate of change of the kinetic energy (E) of the perturbation is governed 

by 

' [v^u'.Sij + R^idu'Jdxj) 2 } dx (1.1) 



where $y — (dUi/dxj + dUj/dxi)/2 is the rate of the strain tensor of the basic flow. 
All quantities are taken to be dimensionless. x = (xx, x%, Xs) = (x,y,z) denotes the 
computation domain. This equation is known as the Reynolds-Orr energy equation. The 
stability problems governed by the energy equation are also referred to non-linear prob- 
lems because the perturbation flow is somehow arbitrary; in particular, its size relative 
to the basic flow is not necessarily infinitesimal. The x-axis or the streamwise direction 
is the direction of the basic flow. We deal with a fluid filling between two solid walls 
situated at ±d w . The origin of the co-ordinates is situated at the mid- way of the walls. 
The basic velocity of the Couette-Poiseuille flow is uni-directional and its normalized 
velocity is given by 

u(z)= { Ml z 1 ) + z, 0^1/2, 

{Z> \ [VA(1-z 2 )+ 2^(1-4)4 1/2<4<1 K ' } 
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for — l^z^l. The Reynolds number R = \J m axd w /v. We assume that u' is spatially 
periodic in the x or y directions so that the integration domain V may extend to inhnity 
in both directions, and the integration can then be taken over exactly one wavelength. If 
the Reynolds number R of a fluid motion satisfies the condition 



R < R n 



mm 



{[f {du'Jdxtfd* I 



}, (1-3) 



then the motion of the basic flow is said to be asymptotically stable in the mean. In the 
above expression, the denominator represents a transfer of energy between the basic flow 
and the perturbation flow. To determine R m in is equivalent to a variational problem. 
The corresponding Euler equation of the variation takes the form of 



R^Au' 



tijSij 



dp'/dxi, dv! i jdxi = 0, u\ — on dV. 



(1.4) 



The last statement is the no-slip boundary condition. The pressure p plays the role of 
the Lagrange multiplier associated with the solenoidal constraint. Let a and j3 stand 
for the disturbance wavenumbcrs in the x and y directions respectively. Consider the 
normal-mode of the form 

K, u' 2 , u' 3 ,p'}(x, y, z) = {/, g, h, q}(z) e iax+i ^, (1.5) 

then the Euler equation reduces to the system of the differential equations 

crLf — U'h = ia q, aLg = i/3 q, aLh — U f = q' , (1.6a,b,c) 
iaf + ifig + ti = 0, /(±1) = g(±l) = h(±l) = (1.6d,e) 

where the operator L = (d 2 /dz 2 — 7 2 ), 7 2 = a 2 + f3 2 and ct" 1 = R/2. The primes stand 
for the differentiation with respect to z. After eliminating the components g and q, we 
obtain 



aL 2 h 



-j 2 U'f-ia{U'h' + U"h), aj 2 Lf = /3 2 U'h + iaa(Lh') (1.7o,&) 



and the boun dary condi t ions f (±1) = h(±l) = h'(±l) = 0. Based on the idea of energy 
dissipation of Reynolds ( 1895b . Orr studied (|1.4|) in an attempt to establish a criterion 
so as to supplement the stability obtained by solving the Orr-Sommerfeld equat ion. H e 



restric ted his calculations to the streamwise disturbances. For plane Couette flow. Uoseph 
( 1966| ) pointed out that the minimum Reynolds number is associated with a spanwise 
distur b ance. Stabili t y crit e ria by the en erg y method have been es tablished in ISharp"el 
(|l905h . lMacCreadiel |l93ll ). lBussel (|l969l ) and I Joseph fc Carmil (|l969h . It should be noted 
that these papers solved the eigenvalue problems whi ch are l inear in a ( or in a 2 ). More 
backg round on the e nergy method can be found in ISerrin (Il959h and IChandrasekharl 
(| 1 9 6 lh . The book by iDrazin fc Reidl (|198lf ) contains a whole section on the role of the 
energy method. In the past three decades, there have been virtually no published works 
on the non-linear stability using the approach of the Reynolds-Orr energy. Our recent 
study shows that there exist some gaps in the knowledge of the stability for the Couette- 
Poiseuille flow. 

The structure of (11.61) is essentially different f rom that of the linearized Navier-Stokes 



equations ( Linl ( 19551 §1) or Drazin fc Reidl (Il981, §25) ). In the linear theory, analysis has 



been facilitated by Squire's theorem, see Sauirei ( 19331 ). However, Squire transform does 
not apply to (11.61) except the case of the linear basic profile. In general we must consider 
fully three dimensional disturbances and consequently we have to deal with some non- 
linear eigenvalue problems. Nevertheless, the mathematical form of (|1.7|) suggests that 
it may be advantageous to convert the equations into the equivalent integral equations. 
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The Nystrom procedures can then be used to efficiently solve the integral equations. Such 
formulations convert the stability problems into the linear or quadratic matrix eigenvalue 
problems. It turns out that use of the integral equations simplifies the analytic treatment 
of the stability problems. In the present paper, we give a summary of the work developed 
from these ideas. 

2. Solution of (fTTTj) for /3 = 

For the two dimensional disturbances in the dire ction of the basic flow, (ll.6[) or (|1.7[) 
is simplified as the differential equation (jOrrl ()l907l )) 

L 2 h(z) = -\(U'ti + U"h/2), h(±l) = h'(±l) = (2.1) 

where A = iaR. The adjoint equation to (|2.1j) takes the form of 

L 2 h\z) = \ {U'h)' + U"tf /2), tf{±l) = h V {±l) = Q. (2.2) 

2.1. An exact solution for plane Couette flow 

We have used two methods to tackle (|2.1j) . The first one is particularly appropriate for 
plane Couette flow. Because of the linearity in the basic velocity profile, the differential 
equation (|2.1j) is first converted into a Volterra integral equation of the second kind. In 
particular, the kernel H of the resulting integral equation is in the convolution form. 
Consequently, the inverse Laplace transform must possess an analytic solution of closed 
form. The four linear independent solutions of the operator L 2 are found to be 

Fi{z) = e az , F 3 (z) = —\{Fi(z)(z + a ) - a F 2 (z)}/2, (2.3a,c) 
F 2 (z) = e~ az , F 4 (z) = -A{a Fx(z) + F 2 (z)(z - a )}/2 (2.36,d) 

where ao = l/(— 2a). Introducing the vorticity rj = Lh, the method of variation of 
parameters is evoked to reduce (|2 . 1 [) to the equivalent Volterra integral equation 

4 fZ 

r,{z)=Y,c j F j (z)+ H(z, S )rj(s)ds (2.4) 

where Cj's are arbitrary constants. The kernel is given by 

H(z, s) = -a a 2 (z - s) ( e a{z - s] - e ~ a{z ^ s] ) . (2.5) 

The special form of the kernel suggests that we seek the solutions of the integral equation 
(|2.4j) by means of the Laplace method. Denote Fj's and H as the Laplace transforms of 
Fj and H respectively. The solution of the integral equation is given by 

where the value of 7 in the Bromwich contour of integration in the complex cr-plane 
is chosen so that all the singularities of the integrand lie on the left of the line, being 
parallel to the imaginary axis. The singularities of the integrand are determined by the 
four roots (assume distinct) of the quartic equation 

er 4 - 2a 2 a 2 + a 4 - Aada = (2.7a) 

where d = —\aa^. Let a^, k — 1, 2, 3, 4 be the roots and they are readily found to be 

0-1,2 = (V6P±Q)/6, ct 3 ,4 = -(VQP±Q)/6 (2.76) 
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where 

P = ^2a 2 + S + 4a 4 JS, Q = J '-(-24 a 2 + 6 S + 24 a 4 /S + 36 dav^/P) (2.7c) 



and 5 = v^a 6 + 27d 2 a 2 + Sadv^a 6 + 81d 2 c?. The inverse transform can be carried 
out by means of partial fraction decomposition and the result is 



/(*) = £<* U(z) + A^> /V^F^ds) 

3=1 V fc=l - 70 / 



(2.8) 



where 6fc = lGofV^ Ylj=i j^ki^k — ^j) 1 - (When R — 4/(3y/3)a, the quartic equation 
admits repeated roots, there will be obvious modifications in the form of the solutions 
and the values of bk) The expression (|2.8j) constitutes an exact solution to the vorticity. 
Once rj{z) is known and so are h(z) and h'{z). The secular determinant, A(a,R) = 0, is 
obtained by applying the four boundary conditions. 

2.2. General velocity profile 
The first two of the linear independent solutions Pj and F 2 remain unchanged while 

F 3A (z) = -\(J Y{z,s)e ±as {aU\s) + U"{s)/2\y s (2.9) 

where Y(z, s) — cto(e az e~ as — e~ az e as ) and F3 takes the upper sign. Then the corre- 
sponding Volterra equation takes an alternative form of 

4 , 

ri{z) = y^c j F j {z) + f (Hi+H 2 )(z,s)Ti(8)ds, (2.10a) 

H 1 (z,s) = - f Z Y( z ,t)^P^U'(t)dt, (2.106) 
Jo ds 

H 2 (z,s) = - I Y(z,t)Y(t,s)U"(t)/2dt. (2.10c) 
Jo 

Clearly the solutions to the integral equation will be much more involved. In the case of 
plane Poiseuille flow, we have 

Hx(z, s) = A{a 2 (z - s)(z + s) - a(z - s) + 1} sinh[a(z - s)]al, (2.11a) 
H 2 (z, s) = {2a{z - s) cosh[a(z -s)]-2 sinh[a(z - s)]}a 2 . (2.116) 

Note that H\ is no longer in the convolution form. 

2.3. Reduction to Fredholm integral equation 

For non-linear velocity profiles, we prefer to use our second method because the self- 
adjoint differential operator L 2 together with its homogeneous boundary conditions has 
a simple inverse. Let the four linearly independent solutions of L 2 be 

01 (2) = e a \ fo(z) = a a {M*)M") ~ Mz)fo(z)} , (2.12o,c) 
fa{z) = e- a \ M?) = «o {M*)M*) ~ M*)M*)} > (2.126,d) 
M?) =z + l, Mz) = Me~ 2a ~ e 2az ), M*) = «o(e- 2QZ - e 2a ). (2.12e,f,g) 
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The Wronskian of the four solutions equals 4a 2 . The Green's function associated with 
L 2 satisfying the homogeneous Dirichlct boundary conditions is found to be 



G(z,y) = a < 



{ {My) + A(o)[^3(i)0 3 (i/) - Mi)Mv)]}M*) - 

{My) + H*)[ik.(i)Mv) - MVMv)]}M*)> 



{Mz) + A(a)[^(l)0 3 («) - *l>iO-)Mz)]}Mv) - 

{M*) + Ha)tyi(i)Mz) - Mi)M*)]}Mv)> y^*, 



(2.13) 



where A(a) = a /{V'2(l)V'3(l) — V'i(l)}- Hence equation (|2.ip is equivalent to a homo- 
geneous Fredholm integral equation of the second kind, namely 



(2.14a) 
(2.146) 



h(z) + X J K(z,y)h(y)dy = 0, 
K(z, y) = U'(y)dG(z, y)/dy + U"(y)G(z, y)/2. 



Note that K(z,y) is non-symmetric. Since G is defined in terms of linear combinations 
of the entire functions, it follows that \K(z, y)\<oo. In addition, the derivative of K with 
respect to y has finite jumps at y = z and therefore it must also be bounded. Accordingly, 
the Fredholm determinant of (|2 . 14[) is an entire function of a and R. Its order is at most 
2/3 for finite a and R. There exists at least one eigenvalue R for every wavenumber a. 
In fact, the spectrum consists of denumerable discrete eigenvalues. 



3. Stability by the linear theory 

Our method of solution can also be applied to the Orr-Sommerfeld equation 

A- 1 L 2 0=([/- c )(0"-a 2 0)-[/"0, 0(±1) = 4>'{±1) = 0. (3.1) 

In the temporal stability problem, c is the complex eigenvalue to be determined via the 
eigenvalue relation A (a, R, c) = 0. The equivalent Fredholm integral equation is 

A~V(*)-/ K L {z,y)<Kv)dy = -c [ K R {z,y)<f>{y)dy, (3.2a) 



J-i J -l 

K L (z,y) = ULG(z,y) + 2U'G'(z,y), K R (z,y) = LG{z,y). (3.26) 

The primes denote the differentiations with respect to y. Briefly, the Orr-Sommerfeld 
equation can be reduced to a generalized matrix eigenvalue problem when the integrals 
are approximated by quadratures. The second derivative of the basic flow does not play 
any role in the formulation. Hence the present method has the advantage in dealing with 
the velocity distributions that only the first derivative is available. It may be useful to 
calculate the critical Reynolds number by the spatial theory (c is taken to be a purely real 
quantity. The wave number a is taken to be complex and is the unknown eigen-modc). 
Both the temporal theory and the spatial theory define the identical marginal stability; 
the critical Reynolds number is same in either theory. In addition, the spatial theory is 
mathematically i ll-pos ed because the velocity and the pressure are unbounded at infinity 



( Drazin fc Reidl (Il98lh l. The ill-posedness applies to most basic flows. The reason why 



one may calculate many spatial modes is that the corresponding temporal modes exist. 
The converse is however not true in general; some eigen-modes of the linear stability are 
well-defined only in terms of the temporal theory. 
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4. Solution of (fTTflt for a = 

For the spanwise disturbances, the differential equation (|1.7p consists in a sixth or- 
der system. The Fredholm integral equation is found by considering the inverses of the 
operators L 2 and L in (|1.7[) . The result is 

a 2 h(z) = -/3 2 J ^ G (z, y) U\y) h(y)dy, (4.1a) 

G (z,y) = J G(z,s)F(s,y)U'(s)ds (4.16) 
and F is the Green's function for the operator L satisfying the boundary conditions. 



5. General solution of (|1.7j) 

The sixth order system of the ordinary differential equations is reduced to the equiva- 
lent Fredholm integral equation 



a 2 h(z) = ~f3 2 / G U'(y)h(y)dy + aia 



l . 



dC i 
G 1 -i 2 G 2 + —U'{y) h(y)dy,(5.1a) 
dy J 



G 1 (z,y) = / d 3 F(z,s)/ds 3 G( S ,y)U'(s)ds, 



-l 
i 



G 2 (z,y)= / dF(z, s )/dsG(s,y)U'(s)di 



(5.1b) 
(5.1c) 



6. Computations 

Practically, the integral equations (|2.14|) and (|4.1|) have been solved by the standard 
Nystrom techniques. We have chosen the iV-point Gauss-Legendre quadrature rule for the 
approximations since the kernels consist of well-behaved smooth functions. Some efforts 
have been given to the dependence of the accuracy of the eigenvalues on the Nystrom 
mesh. We denote the solutions of the integral equations by the epigraph optimization: 

To minimize R = R(a, (3, A) (6.1) 

subject to no constraints. In particular, the solution of (|2.14p or (|4.1[) reduces to the 
minimization problem of a univariate Reynolds-number function for fixed A and to a 
bivariate function for variable A. To avoid excessive computations, we have made use of 
the theory of optimization which avoids any evaluations of the derivative of the Reynolds- 
number function. For every given wavenumber a or /3, the eigenvalue spectrum of R is 
computed. The solutions of (|2.14j) consist of complex conjugate pairs which give identical 
u' 3 while the spectrum of (I4.1[) consists of simple eigenvalues. The variations of R m in 
with the velocity parameter A are plotted in Figure [T] Clearly there is a value of A 
of the basic flow at which R m in is least. The Reynolds-number functions are found to 
be convex and the epigraphs are shown in Figure [2] The main numerical results are 
summarized in Table [TJ In particular, plane Couette flow has been found stable for any 
R < R m in — 44.304 and a min = 1.8934. The symbol, a m i n , denotes the wavenumber 
at which the threshold Reynolds number occurs. The value of R m in agrees with the result 
of Orr (1907) (R — 44.3). Similarly, for purely spanwise disturbances, the numerical 
results are compared in Table [2] The eigen-functions of the spanwise disturbance at 
A = 0.5394 are presented in Figure [3l It is known that the eigenvalue problem defined in 
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Figure 1. Solutions of the Reynolds-Orr energy equation for the Couette-Poiseuille profile (|1.2p . 
Rmin at a m in and at Pmin is shown as a function of the basic flow parameter A. The case on 
the far left-hand side in each plot corresponds to plane Couette flow and at the far right-hand 
side plane Poiseuille flow. 
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Figure 2. Solutions of the Reynolds-Orr energy equations (|2.14|) and (I4.f[) . The least Reynolds 
number R\ of the spectrum at given wave number is shown as a function of a or /3. The dash-dot 
lines are for plane Poiseuille flow, the solid lines for plane Couette flow and the dashed lines for 
the Couette-Poiseuille profile. By the minimization procedure, the lowest point in each curve 
defines the values of Rmin and a m in or /3 m » n . 



Basic Flow ([Of 


A 


Author 


Rmin 




Plane Couette: 





Orr (1907) 


44.3 


1.89 







Present 


44.304 


1.8934 


Plane Poiseuille: 


1.0 


Orr (19071 


87.75 


2.1 




1.0 


MacCreadie 


(1931) 87.6 


2.05 




1.0 


Present 


87.594 


2.0986 


Couette-Poiseuille: 


0.5469 


Present 


42.320 


1.9264 



Table 1. Comparison of the results for streamwise disturbances oc exp(iaa;) 
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Basic Flow JOJ) 


A 


Author 


Rmin 




Plane Couette: 





Pellew & Southwell 


(1940) 20.663 


1.5585 







Present 


20.663 


1.5582 


Plane Poiseuille: 


1.0 


Busse (1969) 


49.604 


2.044 




1.0 


Present 


49.605 


2.0437 


Couette-Poiseuille: 


0.5394 


Present 


19.878 


1.5936 



Table 2. Comparison of the results for spanwise disturbances oc exp(i/3y) 




-1 -0.5 0.5 1 -1 -0.5 0.5 1 -1 -0.5 0.5 1 -1 -0.5 0.5 1 

Figure 3. Solutions of the Reynolds-Orr energy equation (|4.1[) subject to spanwise disturbances 
oc exp(i/3j/). From the left to right, the results of the first four modes of the eigen- functions at 
P min are plotted (.4 = 0.5394). 



(|4.1j) is formally identical to the stability formula tion of a viscous f luid co ntained between 



two rigid boundaries and heated from below, see IChandrasekhar (Il96ll §11). The role of 
the Reynolds number plays that of the R ayleigh number. In additio n, exact eigenvalue 
relations are available (see (46) and (51) of lPellew &; Southwelll ([19401 )). By the procedure 
of the unconstrained minimization, the critical Rayleigh numbers a nd the wavenumbers 
are th en well defined and hav e been calculated. In the scalings of iPellew fc Southwelll 



(|1940T ) or lReid & Harris! (|l958l) . we found the following: 



First even mode: Ra c = 1707.7618, (3 mm = 3.1163, (6.2a) 

First odd mode: Ra c = 17610.3937, p mm = 5.3646. (6.26) 

In comparison to the classical results, good agreement is obtained. The only noticeable 
difference occurs in the wavenumber of the even mode. (The value of f3 m in = 3.117 has 
been quoted in most literature. Although it may appear to be a good numeral, it has 
been obtained by graphical method in the past in an arbitrary way because the R — j3 
graph was very "flat" in the vicinity of the minimum (3.) 

Application of the Nystrom method to (|5.ip leads to a quadratic eigenvalue problem 
of the form 

(<t 2 I + aA + B)h = (6.3) 

where A and B are the discretization matrices, and I the identity matrix and h the 
eigen-functio n vector. The non -linear problem can be solved by means of the linearization 
technique of I Wilkinson (Il965l) . However, it is found that the non-linear problem is not 
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Parameter A R m in ot min p m in 

0.2 20.530 0.0457 1.5627 

0.4 20.150 0.0893 1.5756 

0.6 19.935 0.1315 1.5960 

0.8 22.603 0.1957 1.6417 

1.0 47.196 0.6284 1.8381 

0.5423 19.821 0.1185 1.5889 

Table 3. Results of Rmin for fully three-dimensional disturbances oc exp(iaa; + i/3y) 




completely reducible to the linear problem (|2.14[) as /3— >0 and a being large. Roughly 
speaking, no real Reynolds numbers of (|6.3[) are found for certain values of small /3 and 
large a. Since the results of the calculation for these wavenumbers do not compromise 
the main conclusion, no substantial efforts have been directed to resolve this issue. A 
selection of the computed results is presented in Figure |4j In the plot for Poiseuille flow, 
the minimum Reynolds number is clearly shown to be away from the axis of a = 0, 
indicating the three dimensional nature of the disturbance. The minimization problem 
for the other values of A>0 has been solved and a similar plot to the right-hand plot of 
Figure [1] has been found. The main point of these computations is that every computed 
value of R m in is lower than the corresponding value in the right-hand plot of Figure [TJ 
Some numerical values of the minimum Reynolds number are listed in Table [3l The 
eigen-functions for plane Poiseuille flow at R m in are plotted in Figure [5] The modal 
distributions of the purely two-dimensional and the fully three dimensional disturbances 
are displayed in Figure [5] In every computation, we have examined the first few modes 
of the spectrum because the knowledge of these higher modes enables us to understand 
the critical mode. 



7. Conclusions 

We have shown how the technique of integral equation can be used to solve the 
Reynolds-Orr energy equation. Observation has been made on the mathematical struc- 
ture of the equation for the special case of j3 = 0. It is found that the case of plane Couette 
flow can be solved exactly because the kernel of the integral equation is in the convolu- 
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-1 -0.5 0.5 1 -1 -0.5 0.5 1 -1 -0.5 0.5 1 -1 -0.5 0.5 1 

Figure 5. The eigen- functions of the first 4 modes at a min and fi m i n for plane Poiseuille flow 
subject to fully three dimensional disturbances oc exp(iax + i/3y) . The Reynolds numbers of 
the modes are listed. Subscripts r and i are used to label the real and the imaginary parts 
respectively. 



Streamwise disturbance « e lax Spanwise disturbance « e'^ y 3D disturbance « e" 




1 7 13 19 25 1 7 13 19 25 1 7 13 19 25 
Modes Modes Modes 

Figure 6. The modal structure of the Reynolds-Orr energy equations (|2.14|) . (|4.1[) and (|5.1|) 
for the first 24 modes. In the left-hand side plot, the modal distribution is presented for a m i„. 
The symbols □ are the modes for plane Couette flow, x for the Couette-Poiseuille profile of 
A = 0.5469 and o for plane Poiseuille flow. The similar result is shown in the middle plot for 
fimin, x for A = 0.5394. The results for fully three dimensional disturbances are shown in the 
right plot. The Couette-Poiseuille profile is for A — 0.5423. 



tion form. For general velocity profiles, analytic solutions of th e problem become harder 
to obtain. The minimum Reynolds numbers obtained bv lOrrl ( 1907t ) for plane Couette 
and plane Poiseuille flows have been confirmed and they are numerically correct. For the 
purely two-dimensional disturbances, there exists a particular Couette-Poiseuille profile 
which gives rise to the lowest R m i n for 0^j4^1. For the three-dimensional disturbances, 
the Euler variation problem becomes a non-linear eigenvalue problem for minimizing the 
Reynolds number. It is shown that the minimum Reynolds number is in general smaller 
than the corresponding two-dimensional case for all the Couette-Poiseuille profiles. The 
only exception is plane Couette flow where the minimum Reynolds number is still related 
to the spanwise two-dimensional disturbance in the j/z-plane. For plane Poiseuille flow, 
the minimum Reynolds number of the Euler equation is associated with a fully three- 
dimensional disturbance; R m in is approximately 47.2 with the streamwise wavenumber 



On the Reynolds- Orr energy equation 



11 



being one third of the spanwise wavenumber. The computational results show that the 
Squire's theorem is not valid for the basic flows except the linear profile. 

The p resent study in d icates that certain claims made in the existing literature (for 
instance iDrazin fc Reid ( 1981 )) on the stability characteristics of some basic flows by 
the Reynolds-Orr energy equation are incomplete. 

The threshold Reynolds numbers of the energy criterion are found "rather low" . How- 
ever the criterion has been derived exactly from the Navier-Stokes equations of motion. 
The non-linearity of the equations has been fully captured. On the other hand, the 
Cauchy problem for the Navier-Stokes equations in R 2 (and in all likelihood in R 3 ) has 
been shown to be well-posed. The property of the global regularity has the implication 
that the structural anomalies such as bifurcations and singularities have been ruled out. 
The velocity and the pressure are everywhere unique and smooth during the course of 
the flow development. The solutions of the linear theory (see §2) can be quantitatively 
assessed by direct comparisons with the Navier-Stokes solutions. In fact, some conse- 
quences implied by the linear theory are simply incompatible with the Navier-Stokes 
global regularity. In this note, we are mainly interested in the non-linear stability of the 
energy equation and hence we shall not elaborate on our assessment on the linear theory. 
For the initial data of practical interest, it has been known that the linear theory is 
unable to describe the complete evolution of fluid motion because the linearization fails 
to encapsulate the non-linearity which is the essence of fluid dynamics. 
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